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We perform atomistic computer simulations to study the mechanism of homogeneous dislocation 
nucleation in two dimensional (2D) hexagonal crystalline films during indentation with a circular 
nanoindenter. The nucleation process is governed by the vanishing of the energy associated with a 
single normal mode. This critical mode is largely confined to a single plane of adjacent atoms. For 
fixed film thickness, L, the spatial extent, £, of the critical mode grows with indenter radius, R. For 
fixed R/L, £ grows roughly as £ ~ L 0A . We, furthermore, perform a mesoscale analysis to determine 
the lowest energy normal mode for mesoscale regions of varying radius, r mcso , centered on the critical 
mode's core. The energy, A mcso , of the lowest normal mode in the meso-region decays very rapidly 
with r meS o and A me so ~ for r mcso > £. The lowest normal mode shows a spatial extent, £meso, 
which has a sublinear power-law increase with r moso for r mcso < £ and saturates at r moso ~ 1.5 £. 
We demonstrate that the £ m cso/£ versus r meS o/£ curve is universal, independent of film thickness 
or indenter radius. The scenario that emerges is one where the analysis of small regions, r mcso < £, 
in the material can reveal the presence of incipient instability even when the region being probed 
is much smaller than the spatial extent of the critical mode. However, the mesoscale analysis 
gives good estimates for the energy and spatial extent of the critical mode only for r mcso > 1.5 £. 
In this sense homogeneous dislocation nucleation should be understood as a quasi-local phenomenon. 

PACS codes: 61.72.Lk, 62.20.-x, 62.25.-g, 81.40.Lm. 



I. INTRODUCTION 

It has long been understood that plasticity in single 
crystals is governed primarily by the motion and produc- 
tion of dislocations 1 . Typical crystals contain many 
pre-existing defects that can act as sources for the pro- 
duction of new dislocations. However, in nanoscale sam- 
ples, one can obtain essentially defect free crystals. When 
subjected to inhomogeneous loading, for example under- 
neath a nanoindenter, new dislocations will be created 
in these perfect crystals. The question of where in the 
sample and under what conditions this happens is still 
surprisingly contentious. 

One naive expectation is that the production of new 
dislocations on a particular slip system is governed by 
the resolved shear stress, much in the same way that 
the motion of a pre-existing dislocation is known to be 
governed by the resolved stress [2] . Li and co-workers [3J 
first pointed out that this simple and intuitive idea is not 
correct and demonstrated this quantitatively in atom- 
istic models. Instead, they proposed to use a modified 
version of a stability criterion first given by Hill in the 
context of a continuous elastic medium [3], the so-called 
A-criterion [5]. 

Later, Miller and co-workers [H [7] showed that, al- 
though the A-criterion worked for some potentials and 
some indentation geometries, it failed for others. The 
A value could indicate instability at many locations in 
the material even while the system remained mechani- 
cally stable. Furthermore, when the system did become 
unstable, the motion responsible for triggering the insta- 
bility, the so-called dislocation embryo, was localized in 
a region far from those atoms considered least stable in 
terms of the A-criterion. Miller and Acharya [B] initially 



proposed a new stability criterion derived from Field Dis- 
location Mechanics (FDM) theory, but later Miller and 
Rodney (MR) {7\ showed that this criterion had its own 
shortcomings. 

Most recently, MR [7] proposed to examine the en- 
ergy eigenmodes of regions of space extended over many 
atoms, but smaller in spatial extent than the full system. 
In particular, since the dislocation embryo was found to 
be localized largely on a disk-like region, they chose disk- 
like regions of adjacent planes of atoms a posterior to 
have a large enough spatial extent to encompass a whole 
embryo. They found that the lowest energy eigenvalue 
of the disk-like region was an excellent predictor of the 
location of the embryo. Since then, others have used 
similar approaches to study various types of nucleation 
phenomena [5HTU], 

This initial work left several important issues open. 
On one hand, MR showed evidence for the growth of the 
embryo with the radius of curvature of the indenter tip. 
On the other hand, they did not discuss how the choice 
of size of the disk-like meso-region would affect the anal- 
ysis. One might wonder, for instance, if a meso-region 
would be able to detect the instability if it were much 
much smaller than the intrinsic embryo. An important 
related question is how rapidly the lowest energy eigen- 
value approaches zero with increasing meso-region size. 

In this work, we show, in agreement with MR, that 
the intrinsic dislocation embryo grows with the radius 
of the indenter. However, here, we focus on films of fi- 
nite thickness and show that even when the ratio of film 
thickness to indenter radius is kept constant, the embryo, 
surprisingly, grows with the film thickness. One can ar- 
gue dimensionally that this effect would be completely 
absent in any continuum description, as in reference [3J, 
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that does not contain an additional length scale in addi- 
tion to the film thickness and indenter radius. We also 
show that a meso-scale analysis can detect the presence 
of a dislocation embryo even when the meso-region is 
much much smaller than the embryo itself. In particular, 
we show that the inferred spatial extent of the embryo 
captured by the meso-region becomes precisely equal to 
the true spatial extent of the embryo, but only after the 
meso-region fully encompasses the embryo. 

The rest of this paper is organized as follows. In sec- 
tion [Tl] we describe the details of our simulation setup 
and our athermal, quasitatic indentation protocol. In 
section |III| we give a give a brief kinematic description 
of the mechanism of homogeneous dislocation nucleation. 
In addition, we precisely define our measure for the em- 
bryo size, £, and present our results on its scaling proper- 
ties. Section HVl contains our rationale and results on the 
mesoscale analysis. In particular we discuss how the low- 
est meso-scale eignevalue and the inferred spatial extent 
of the dislocation embryo scale with meso-region size. Fi- 
nally in section [V] we conclude with a brief summary and 
discussion on implications of our results and an outline 
of future work. 



II. SIMULATION PROTOCOL AND 
FORMALISM 

We perform nano-indentation simulations via energy 
minimization dynamics. The LAMMPS molecular dy- 
namics framework [11) is used to perform the simulations. 
The system is two dimensional (2D) with a Lennard- 
Jones (LJ) interaction potential. All energies and dis- 
tances are measured in Lennard- Jones units. We use 
a cutoff of 2.5 for the interactions. Minimization stops 
when no force on any particle exceeds 10 -8 . 

The system is periodically replicated on the sides and 
bounded by a half-space of immobile particles at the bot- 
tom that interact with bulk particles through identical LJ 
interactions. The nano-indenter is modeled by a circular 
arc of L J particles (with spacing equal to the lattice con- 
stant) that are constrained to move as a rigid body. The 
hexagonal crystal is prepared with first nearest neighbors 
placed at the equilibrium bond spacing and subsequently 
allowed to relax. The initial quenched stresses are much 
smaller than the stresses induced by the subsequent in- 
dentation. The three crystal axes are placed at ±n7r/3 
with respect to horizontal. 

For our 2D system there are three geometrical param- 
eters: the crystalline film's width and thickness and the 
indenter radius. For the results in this paper, the sys- 
tems are all of aspect ratio (ratio of width to thickness) 
1:1. We verified that if wc double the width while keep- 
ing the same film thickness (an aspect ratio of 2:1), our 
results don't change. We denote the film thickness by L. 
In figure [T] we sketch the three geometries studied: for 
each L, we perform simulations for 3 different indenter 
radii, R, which are equal to L, 2L and 3L. 



FIG. 1: Schematic of our 3 indenter-film geometries. The 
width and depth in each case are equal to L. Left: iiindenter = 

L. Center: indenter = 2L. Right: iiindenter = 3L. 




FIG. 2: Top: Potential energy, U, of the crystalline film as 
function of indenter depth. Bottom: Load, F, on the indenter 
atoms in the vertical direction as a function of indenter depth. 



Following standard procedures, we simulate the ather- 
mal, quasi-static (AQS) deformation by alternately mov- 
ing the indenter by a discrete step and minimizing the 
energy [IS]. When the system loses mechanical stability 
and a dislocation loop (or edge dipole in 2D) is nucle- 
ated, the energy and indenter load will decrease at fixed 
indenter depth. This energy decrease defines an event. In 
figure [2] we show the total potential energy of the system 
and the load on the indenter during a typical AQS sim- 
ulation as a function of the indenter depth, D. At the 
depth, D = 0, the indenter tip just reaches the undis- 
turbed top surface of the crystal. The system loads and 
gains energy reversibly during indentation until a drop in 
the energy at constant indenter depth that corresponds 
to an irreversible dislocation nucleation event. To study 
the nucleation process in detail, we restart the simulation 
in a configuration just before the nucleation event and 
restart the stepping protocol with decreased step size to 
converge to the nucleation point with a precision limited 
only by our energy minimization algorithm. 

The second derivative of the energy with respect 
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FIG. 3: Smallest 4 energy eigenvalues, A p , of one of our 
systems (L = 45 and .Rindenter = 3L), as function of 8D = 
D c — D. The dashed (cyan) line tracks the critical normal 
mode along which the system is driven to instability. 

to particle positions, or Hessian matrix, is defined as: 

n2ir 

Hiajp = dx g x ■ We use latin to index particle number 
and greek to index cartesian components. At any partic- 
ular indenter depth we compute the few smallest energy 
eigenvalues of the Hessian (and corresponding eigenvec- 
tors) using Krylov subspace based routines in the MAT- 
LAB eigsQ function with default parameters. 



III. KINEMATIC DESCRIPTION OF 
DISLOCATION NULCEATION 

Energy minimization dynamics cause the system to fol- 
low local minima of the atomistic potential energy surface 
(PES) as they move in configuration space under changes 
in D. The PES itself changes smoothly with D, however, 
topological changes in its stationary points occur. In par- 
ticular, incrementing D may cause flattening of the PES 
along some direction |13j . In the simplest scenario, the 
system suffers a mechanical instability, such as a homoge- 
nous dislocation nucleation, when the PES has a single 
normal mode with a vanishing eigenvalue. Generically, 
to lowest non-trivial order, the potential energy surface 
(PES) of an atomistic system near the loss of mechanical 
stability should be given by a cubic term in the reaction 
co-ordinate plus a term which is bi-linear in the reaction 
co-ordinate and the external control parameter, D [14] . 
Such a scenario is known as a saddle-node bifurcation 
or fold catastrophe in control theory. The bifurcation 
should be of saddle-node type as long as there are no 
symmetries in the system and there is only a single de- 
gree of freedom in the high dimensional landscape along 
which the curvature of the PES vanishes as the system 
is driven externally. 

In our simulations we see that the system is driven to 
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FIG. 4: Top: Normal mode corresponding to lowest eigen- 
value at: SD = D — D c ~ 10 -6 , for the system L — 45, 
R — 3L. Center: £1 field computed from the normal mode 
with i w 27r/3. 



instability along a single normal mode and the loss of 
mechanical stability coincides with nucleation of a sin- 
gle pair of dislocations [TS]- We denote D c , the critical 
depth, as the depth of the indenter where a dislocation 
pair nucleates. In figure [3j we show the lowest energy 
eigenvalues of one of our systems as a function of SD 
(where SD = D c — D). As expected from saddle- node 
bifurcation theo ry o ne eigenvalue descends through the 
spectrum like yj SD. This behavior is typical of all our 
systems. 

We call the normal mode corresponding to the eigen- 
value vanishing as \J SD, the critical mode. At the bi- 
furcation point, this critical mode points, in the configu- 
ration space, along the reaction coordinate along which 
the system is driven to instability. Although saddle-node 
bifurcation theory gives no indication of the spatial struc- 
ture of the critical mode, it tells us when the lowest eigen- 
mode in the system gives an accurate representation of 
the reaction co-ordinate for the pair nucleation event, 
i.e., the dislocation embryo. For the system under study 
in figure |3J the direction on the PES that is flattening - 



corresponding to the embryo - only corresponds to the 
lowest eigenmode for SD ^ 10 -3 and the lowest normal 
mode is only a meaningful proxy for the embryo below 
this SD. In the analyses and results below, we always 
verify explicitly that we are close enough to the bifur- 
cation point so that the lowest eigenvalue is descending 
through the spectrum according to the saddle-node bi- 
furcation analysis. 

Figure [4] (top), shows the critical mode in a typical 
small system. Dislocation dipole nucleation is charac- 
terized by the anti-parallel motion of a small number of 
atoms on adjacent crystal planes. As reported by MR, 
we find that the embryo involves more than the neigh- 
borhood of a single atom: it is nonlocal. Following the 
spirit of MR but departing in the particular details, we 
quantify the jump in the mode vector by the following 
procedure. First we triangulate the lattice. Then, on 
each triangle, we make a linear interpolation of the mode 
vector defined at each of the three nodes. Next we select 
a particular crystal axis to analyze and choose only the 
component of the mode vector along that axis. Finally 
we take the derivative of that quantity in the direction 
normal to the crystal axis of interest. This quantity is 
called SI. 

Figure [4](top) shows a typical critical mode, and fig- 
ure [^bottom), shows the associated Cl field correspond- 
ing to the 8 = —it/3 crystal axis. As expected from fig- 
ure Qtop), f2 is practically zero everywhere except in the 
core of the embryo. The embryo center is defined as the 
centroid of the triangle with the greatest value of fl. We 
define the slip-line as the line along the axis of interest 
containing the embryo center. Given the extended region 
of non-zero fl along the slip-line, it does not make sense 
to talk of any point of nucleation. Rather, the nucleation 
event is spatially extended. 

In figure [5] we extract and plot fl(s) profiles along the 
slip-line, for systems of thickness L = 45, 90, 135, 180, 
225, 270 and for all 3 indenter-crystal geometries (shown 
in figure [lj. s indicates distance along the slip- line with 
s = indicating the point of maximum CI. s > corre- 
sponds to locations closer to the surface and s < cor- 
responds to locations closer to the substrate. The over- 
all normalization of the critical eigenmode is arbitrary, 
and the f2 curves are all normalized by their maximum 
value. Figure [5] (top) gives fl(s) for various overall film 
thickness, L, but for a fixed indenter-crystal geometry 
(R = 2L). Figure [5] (bottom) gives Sl(s) for various in- 
denter radius at a fixed film thickness of L — 180. 

The fl(s) profiles of all systems show smooth decay 
from their maximum. For all L and R they are well 
approximated by a Gaussian profile near the center. To 
define an embryo size, we fit the f2(s) profiles, from their 
peak to half their maximum value, to Gaussian functions 
of the form: 
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It is apparent from the Q(s) profiles of figure [5] that £ 
increases with L for fixed R/L and £ increases with R at 



FIG. 5: Q(s) curves on slip-planes, scaled by their maximum 
value f2 max . Top: For different different film thicknesses L 
and one geometry, -Rmdenter = 2L. Bottom: For L — 180 and 

R indenter — 



fixed L. 

Before looking at the variation of £ with L, it is instruc- 
tive to look at the scalings of a couple of more elementary 
quantities. We denote the depth of the defect-center (the 
centroid of the nearest neighbor triad with the highest f2) 
measured from the top surface of the undisturbed crystal 
as Y* . The indenter depth of the last stable configuration 
measured from the same reference as before is denoted by 
D*. Figure [6] shows £>* and Y* scaled by film thickness, 
L, for all three indenter-crystal geometries. To a very 
good approximation these quantities scale linearly with 
L. This simple linear scaling seems to suggest that a 
length-scale free continuum model might suffice to de- 
scribe dislocation nucleation. However we show below 
that £ scales in a much more complicated way with L 
than and Y*. 

In figure [7| we present the value of £, scaled by film 
thickness, L, for each of our 3 indenter-crystal geome- 
tries. We estimated £ using the generalized linear least 
squares method [TB], to log Cl(s) curves upto half their 



5 



0.1- 



% 0.08 - 



0.06 



0.4- 



0.3- 



0.2- 



50 



^ i 1 r 
-© o 



o^)R= 


=L 




=2L 


R= 


=3L 



e 

j , L 



G-O R= 


=L 


mR= 


=2L 


R= 


--3L 



100 



150 200 250 300 
L 



FIG. 6: Top: Depth of indenter in the last stable config- 
uration, D*, scaled by film thickness, L, for all geometries.. 
Bottom: Depth of the defect center, Y*, scaled by film thick- 
ness, L, for all geometries. (All depths are measured from the 
undeformed top surface of the crystal.) 



IV. MESOSCALE ANALYSIS OF INCIPIENT 
DISLOCATION 

A. Rationale for mesoscale analysis 

A system about to nucleate a dislocation pair has a 
critical normal mode which has strong jumps (high fl) 
in only a small portion of the crystal. This motivates 
mesoscale analysis of the incipient defect where we 'freeze 
out' those degrees of freedom outside some some circu- 
lar meso region and analyze its linear elastic properties. 
In particular, we compute the lowest energy eigenvalue, 
A m j n , and eigenmode of meso-regions of increasing radii 
centered at the embryo, in configurations close to insta- 
bility (SD ~ 10 -6 ). MR showed that the lowest eigen- 
value of a meso-region large enough to completely encom- 
pass the true embryo is an excellent predictor of homoge- 
neous dislocation nucleation [7j. Here, we systematically 
investigate what information can be gleaned from analy- 
sis of meso-regions smaller than £. 
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FIG. 7: £/L as function of film thickness, L, for all 3 
indenter-crystal geometries obtained from Q(s) curves. 
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maximum value. We have checked that our results were 
not sensitive to the details of the fitting procedure. To a 
good approximation, our results show that for fixed R/L, 
£/L ~ L~ os . The prefactor of these curves increases 
with indenter radius, as has been previously shown in 
the limit of L S> R 7] ■ The nonlinear scaling of £ with L 
implies that no length-scale free continuum model, sug- 
gested by the simple linear scaling of and D*, can 
possibly predict the spatial structure of the embryo. 



FIG. 8: Top: Characteristic A m i n versus r mcso curves for 
an undeformed and a deformed crystal close to nucleation 
(8D ~ 10~ 6 ). Bottom: Meso region normal modes corre- 
sponding to the lowest energy eigenvalue of meso-region with 
radius r mcso = 8: from an undisturbed crystal (left) and a 
configuration close to dislocation nucleation (right). 

For a homogeneous, linear, elastic continuum with 
fixed boundaries, the energy eigenmodes scale as 1/L 2 
(where L is a characteristic measure of system length). 
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figure |3j the plateau height will go to zero as v SD. The 
crossover r meso increases with L in accord with figure [7] 
It may be tempting to think of the exponential decay of 
Amin with r mcso as intrinsic property 9 and to extract 
a length scale associated with the embryo. However, it 
should be kept in mind that the plateau value of A m ; n 
is a function of SD and not an intrinsic property of the 
system. Systems prepared closer to the bifurcation point, 
at smaller SD, have a lower plateau value and a more 
rapid decay to the plateau, so there is nothing intrinsic 
about the decay rate. 



FIG. 9: ln(A m i n ) versus r mcso curves for various film thick- 
nesses, L, and one of our geometries with -Rmdcntcr = 2L. 



This is true for atoms in a crystal as long as the domain 
size is large enough that continuum effects dominate the 
behavior of long wavelength normal modes and the crys- 
tal atoms are confined to small displacements. In the top 
panel of figure [8] we demonstrate the validity of this ap- 
proximation for an undeformed crystal by showing that: 
~ iAmcso- However, the same plot shows that, 
for a typical deformed configuration close to nucleation 
(SD ~ 10~ 6 ), A m ; n falls off faster than any power law 
with increasing meso-region radii. |19] 

For a given meso-region of radius r meso , the normal 
modes corresponding to A m ; n for the undeformed config- 
uration and the deformed configuration about to nucleate 
a dislocation are very different. The bottom image of fig- 
ure |8j shows that the lowest energy normal mode for the 
undeformed configuration, at r mcso = 8, is a long wave- 
length mode with dominant wavelength equal to roughly 
twice the diameter of the meso-region (as expected from a 
continuum approximation). However for the same r m eso! 
the adjacent image shows that the lowest normal mode 
for a typical deformed configuration has a Gaussian em- 
bryo resembling the critical normal mode shown in fig- 
ure|4](top). This suggests that the structure of the lowest 
meso-region normal mode has information about the true 
dislocation embryo, even when the eigenvalue is far from 
zero. 

Once r mcso becomes large enough that the entire em- 
bryo is encompassed, one gains no new information by 
considering larger meso-regions. Hence we expect that 
the A m i n versus r mGSO curve for a system near nucleation 
must plateau at some r mGSO . The crossover to the plateau 
must be dependent on the embryo size £. In figure [9] we 
show ln(A m j n ) versus r mcso curves for one of our represen- 
tative indenter-crystal geometries (R = 2L), for configu- 
rations close to nucleation (SD ~ 10 -6 ). We see that for 
deformed configurations close to nucleation, A TO j„ rapidly 
decays with r meso before plateauing. In each case, there 
is a roughly exponential drop of A m ; n before it plateaus 
at some large r meS o- 

The plateau height, the value of A min in the limit of 
large r meS o, depends only on SD: as we have shown in 



B. Mesoscale analysis of critical mode structure 

Given that dislocation nucleation is a nonlocal phe- 
nomenon with a spatially extended embryo, the central 
question we address with mesoscale analysis in this pa- 
per is: how big a region in a crystal is big enough for 
detecting an impending nucleation event and the asso- 
ciated embryo? Unlike earlier work [345] that focussed 
on the lowest energy eigenvalue, A m i n , of meso-regions, 
we analyze the the spatial structure of the lowest nor- 
mal mode of the meso-region. The main reason for our 
choice is that the height of the plateau of A m ; n versus 
''meso depends very sensitively on SD whereas the spatial 
structure of the critical normal mode does not. Again, in 
all our analysis we have been careful to choose configura- 
tions where SD is small enough that the lowest eigenmode 
of the full system corresponds to the dislocation embryo. 

First, in figure [lO] we show the lowest normal mode 
in a meso-region for a typical system near its nucleation 
threshold. For comparison, the true critical mode for 
the full system is also shown. We indicate the usual 
unconstrained eigenvalue analysis of the full system by 
r mcso = oo. The lowest normal mode in each meso-region 
has a Gaussian fi(s) profile like the true critical mode. 



Note that for the system studied in figure 10 



is smaller than the spatial extent of the slip-plane £ ~ 13 
measured from the £l(s) profile of its critical mode. So, 
we conclude that even meso-regions much smaller than 
the intrinsic dislocation embryo capture the center and 
essential structure of the true critical mode. 



In figure 11 we show f2(s) profiles for the lowest modes 
of meso-regions of increasing size. These curves converge 
toward the Q(s) profile of the true embryo as r moso grows. 
To quantify the spatial extent of the lowest mode in the 
meso-regions, in precisely the same way as for the true 
critical mode, we fit each fl(s) profile to a Gaussian to 
extract a length scale, £ meS o- We label the embryo size 
obtained from the true critical mode in the full-system 
eigenmode analysis as £oo. 

Figure [12] shows £ meS o versus r mcso for a representative 
indenter-crystal geometry (R — IV). As we might ex- 
pect from figure [ll] for small r meS o > £meso grows before 



plateauing at £00. The increase of £ meS o with 



for 



meso regions smaller than ^ x approximately follows a 
power law. This increase is expected in the small r mcso 
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FIG. 10: Structure of the lowest energy eigenvalue normal 
mode for 3 different r moso from the system with L = 135, 
R — 3L and SD ~ 10~ 6 . r mcso = oo corresponds to the entire 
system. Note how quickly the meso normal mode captures 
the structure of the incipient defect. 



FIG. 12: Top: £ mcso as a function of r mcso for various film 
thicknesses, L, and one of our geometries with -Rmdentcr = 
2L. Inset shows that the curves can be made to collapse by 
rescaling the axes with their plateau value £oo. 




FIG. 11: Q(s) curves on slip-plane computed from lowest 
normal modes of various mesoscale regions from the same 
system as figure [l0| 



regime where r meS o < £oo- Recall that £oo grows with L 
as shown in figure [7] 

In the inset of figure [12] we show the rescaled curves 
6ncso/£oo versus fmcso/^oo. The collapse of the rescaled 
curves is also excellent for the other indentation geome- 
tries, R = L and R = ZL (not shown). Furthermore, in 
figure [13} we show three rescaled curves, one from each of 
the three indentation geometries at L — 150. The ^ mcso 
vs r meso curves, upon rescaling by £oo, is universal and 
does not depend on R or L. The solid black line corre- 
sponds to £ m0 so/£oo — 1 an d the dashed black line is the 
£meso = r mcso guide line. The increase of £ n 
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FIG. 13: Collapsed £ meS o versus r meS o curves (rescaled by 
^oo) for systems with L = 180 but the 3 different indenter 
geometries (R = L, R = 2L and R = 3L), showing the ex- 
istence of a universal curve. £ mcso initially grows sublinearly 
with r moso before plateauing at £oo 



is slower than linear for small meso-regions. So, although 
even small meso-regions host a lowest mode that resem- 
bles a dislocation embryo, is only for r mcso > 1.5£oo, that 
the meso-region gives a reliable characterization of the 
embryo size. It is in this sense that we say that nucle- 
ation is a quasi-local phenomenon. 



V. DISCUSSION AND SUMMARY 

In this paper we studied the nucleation of dislocation 
dipoles in the bulk of perfect 2D crystals subjected to 
nanoindentation with a circular, atomistic indenter under 
athermal, quasistatic conditions. We first described the 
kinematics of homogeneous dislocation nucleation, show- 
ing that nucleation involves driving a system to instabil- 
ity along a single, critical, normal mode. We demostrated 
the presence of scalings expected of such a saddle-node 
bifurcation type instability. The critical normal mode 
was shown to be nonlocal with strong jumps, largely con- 
fined to a pair of adjacent crystal planes. We introduced 
a measure, £, for the spatial extent of the embryo and 
showed that for fixed indenter radius, i?, £ grows with 
L, while for fixed R/L, £ ~ L 0A . This universal scaling 
has an important consequence. On simple dimensional 
grounds, no continuum theory containing only R and L 
as geometrical parameters (as used in, e.g., reference [3] 
can be used, even in principle, to describe the embryonic 
size scaling. 

We then performed a mesoscale analysis of configura- 
tions at the stability threshold, showing that small, non- 
local regions of the crystal, centered at the embryo core, 
contain significant information about an incipient nucle- 
ation event. However, unlike previous work that utilized 
the minimum eigenvalue of the meso-regions as the main 
analysis tool [THS], we focused our attention on the spa- 
tial structure of the lowest meso-region normal mode. We 
found that the relation between £ mcS o, the embryonic size 
inferred from the meso-region, and r meso , the size of the 
meso-region itself is universal. The lowest normal mode 
and eigenvalue were found to provide excellent estimates 
of the structure and energy of the true critical mode, but 
only for meso-regions larger than r mcso > 1.5 £. This 
scenario leads us to think of homogeneous dislocation 
nucleation as quasi-local: full information about the na- 
ture of the embryo can only be obtained by analyzing 
sufficiently large regions, however, its existence can be 
inferred by examining regions much much smaller than 



its intrinsic size. 

The ultimate use of the analysis presented here would 
be to inform coarser-grained models, that do not explic- 
itly take into account the atomic degrees of freedom, 
about the creation of new dislocations out of the void. 
For example, in field dislocation mechanics (TTJ [TH] one 
introduces a continuous field to represent the dislocation 
density. It is our hope that a criterion for dislocation 
nucleation based on a meso-scale analysis like we pre- 
sented here could serve as a guide for the introduction of 
atomistic details at the dislocation embryo in concurrent 
multi-scale schemes built on field theories like FDM. 

In the future we will check to see if the basic picture put 
forward here remains the same in other crystallographic 
orientations and in 3D with more realistic interaction po- 
tentials for various metals. In the spirit of building a 
practical criterion for predicting homogeneous disloca- 
tion nucleation, an obvious extension of our mesoscale 
analysis is to consider meso-regions centered at different 
points in the crytsal and not just at the embryo center, 
since it is not known a priori where a dislocation will 
nucleate. Our preliminary results (not presented here) 
indicate that a meso-region need only have some overlap 
with the true embryo core in order to have a lowest nor- 
mal mode indicate the presence of the embryo and the 
location of its core. Although our previous results [T5] 
showed that the approach to nucleation does not apprea- 
ciably affect the structure of the critical normal mode, 
more work is needed to understand the effects of the 
proximity to nucleation, SD, on mesoscale analysis, and 
an initial study of the impact of 8D is already underway. 
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